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I.  INTRODUCTION 


A  liquid-filled  projectile  can  become  unstable  in  flight  when  resonance 
occurs  between  the  angular  motion  (nutation)  of  the  shell  and  certain  non- 
axisymmetric  inertial  oscillations  of  the  spinning  liquid.  Theoretical  deter¬ 
minations  2’ 3  of  frequencies  of  these  oscillations  (eigenfrequencies),  to¬ 
gether  with  their  associated  decay  rates,  have  been  made,  mainly  for  liquids 
in  solid  body  rotation  after  completion  of  spin-up.  These  stem  from  the  work 
of  Stewartson 4  and  Wedemeyer.5 > 6  Reference  7  describes  flow  field  pressure 
measurements  made  to  determine  eigenfrequencies  experimentally  for  solid  body 
rotation,  and  Reference  8  treats  measurements  made  during  spin-up  from  rest. 

The  next  step  is  to  determine  the  liquid  moment  acting  on  the  casing  and 
then  to  predict  the  angular  motion  of  the  projectile.  Because  of  a  simplify¬ 
ing  approximation*  the  early  predictions  of  pressure  moment 2  were  limited  to 


*  The  sum  of  a  Laurent  series  is  replaced  by  a  single  term;  see,  e.g.,  Eq. 
(5.10)  of  Reference_  4. 

1. '  J.  T.  Frasier  and  W.  E.  Scott,  "Dynamics  of  a  Liquid-Filled  Shell,  "  BRL 

Report  No.  1391,  February  1968.  AD  667365. 

2.  Engineering  Design  Handbook,  Liquid-Filled  Projectile  Design,  AMC  Pamphlet 
706-165,  April  1969.  AD  853719. 

3.  C.  W.  Kitchens,  Jr.,  N.  Gerber,  and  R.  Sedney,  "Oscillations  of  a  Liquid 
in  a  Rotating  Cylinder:  Part  I.  Solid-Body  Rotation,  "  BRL  Technical 

'  Report  ARBRL-TR-02081,  June  1978.  AD  A057759. 

4.  K.  Stewartson,  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid,  " 

J.  Fluid  Mech.,  Vol.  5,  Part  4,  September  1959,  pp.  577-592. 

5.  E.  H.  Wedemeyer,  " Dynamics  of  Liquid-Filled  Shell:  Theory  of  Viscous 
Corrections  to  Stewartson' s  Stability  Problem,  "  BRL  Report  1287,  June 
1965.  AD  472474. 

6.  E.  H.  Wedemeyer,  "Viscous  Corrections  to  Stewartsons1  s  Stability 
Criterion,  "  BRL  Report  No.  1325,  June  1966.  AD  489587. 

7.  R.  D.  Whiting,  "An  Experimental  Study  of  Forced  Asymmetric  Oscillations  in 
a  Rotating  Liquid-Filled  Cylinder,  "  BRL  Technical  Report  ARBRL-TR-02376, 
October  1981.  ADA  107948. 

8.  S.  Stergiopoulos,  "An  Experimental  Study  of  Inertial  Waves  in  a  Fluid 
Contained  in  a  Rotating  Cylindrical  Cavity  During  Spin-Up  From  Rest,  " 

PhD.  Thesis,  York  Univeristy,  Toronto,  Ontario,  February  1982. 
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nutational  frequencies  lying  close  to  liquid  eigenfrequencies.  Two  recent  "  f 

studies9*10  have  produced  liquid  moment  calculations  valid  for  all  angular  ' 

frequencies.  Reference  9  treats  only  the  pressure  moment  on  completely-filled  § 

projectiles  using  viscous  perturbation  equations.  Reference  10  treats  both  f 

presssure  and  viscous  shear  moments  for  partially  and  totally  filled  projec-  f 
tiles,  with  and  without  central  rod;  this  work  is  an  extension  of  the  J 

Stewarts on-Wedemeyer  theory  which  employs  the  inviscid  perturbation  equations.  f 
In  this  report  we  extend  the  analysis  of  Reference  9,  which  uses  the  viscous  | 

perturbation  equations,  to  include  the  shear  moment;  again,  we  treat  only  the  i. 

filled  shell .  | 

There  are  four  basic  assumptions  in  all  the  studies,  including  the  present  | 
one:  (1)  The  angle  of  yaw  is  very  small,  permitting  linearization  of  Navier-  * 

Stokes  equations  and  boundary  conditions.  (2)  The  projectile  is  traveling  in  1 

a  straight  trajectory,  is  nutating  at  a  constant  rate  about  a  point  on  its 
axis,  and  experiencing  exponential  yaw  growth  with  time.  (3)  The  initial 
state  of  the  liquid  is  solid  body  rotation  at  a  spin  rate  that  remains  un¬ 
changed  even  after  the  perturbation  is  applied.  (4)  The  timewise  variation 
of  the  flow  variables  is  the  same  as  that  of  the  motion  of  the  shell.  ? 

Gyroscope  experiments11,  12  have  provided  simultaneous  measurements  of  | 

(1)  nutational  frequency*  t  and  (2)  yaw  growth  rate  e,  the  two  param-  | 

eters  that  describe  the  angular  motion.  Theoretical  outputs  will  be  compared  1 

with  these  results.  * 


II.  ANGULAR  MOTION  OF  PROJECTILE  ' 

Here  we  summarize  Chapters  II  and  VI  of  Reference  9.  Two  coordinate  sys-. 
terns  are  considered.  The  first  is  an  inertial,  earth-fixed  system  of  axes  x, 
y,  z.  The  x-axis  coincides  with  the  projectile  velocity  vector,  and  the  z- 
axis  lies  in  the  vertical  plane;  then  the  y-axis  is  directed  so  as  to  form  a 
right-handed  system.  The  second  system  is  the  aeroballistic  x,  y,  z  non- 


*  Definitions  of  terms  are  given  in  LIST  OF  SYMBOLS  Section.  ■ 

9.  N.  Gerber ,  R.  Sedney,  and  J.M.  Bartos ,  "Pressure  Moment  on  a  Liquid-Filled  : 
Projectile:  Solid  Body  Rotation,"  ARBRL-TR- 02 422,  October  1982.  (ADA  120567) .\ 


10.  C.H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  With  a  Viscous  £ 

Liquid  Payload,"  BRL  Memorandum  Report  ARBRL -MR- 3194,  July  1982.  t 

AD  A118676.  '  '  i 

| 

21.  W.P.  D’Amico,  Jr.,  and  T.H.  Rogers,  "Yaw  Instabilities  Produced  by  £ 

Rapidly  Rotating,  Highly  Viscous  Liquids,"  AIAA  Paper  81-0224,  AIAA  19th  | 

Aerospace  Sciences  Meeting,  St.  Louis,  Missouri,  12-15  January  1981.  ? 

12.  R.  Whiting  and  N.  Gerber,  "Dynamics  of  a  Liquid-Filled  Gyroscope:  Update  ; 

of  Theory  and  Experiment,”  BRL  Technical  Report  ARBRL-TR-02221 ,  March 

1980.  AD  A083886.  f 

1? 
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rolling  system  which  has  the  x-axis  il ong  the  projectile  axis  of  symmetry,  and 
the  z-axis  initiallyjn  the  vertical  plane.  These  systems  are  shown  in  ’• 
Figure  1;  the  y  and  z  axes  are  omitted  for  clarity. 

The  x  =  0  and  x  =  0  values  are  located  at  the  midplanes  of  the  unyawed  and 

yawed  cylinders,  respectively.  The  x-axis  is  nutating  about  the  x-axis  at  the 
angle  =  K^(t).  The  components  of  the  projection  in  the  y,  z  plane  of  a 

unit  vector  lying  on  the  x-axis  are  denoted  by  ny£  and  n^. 

The  yawing  motion  is  characterized  by  two  variables,  a  and  ft.  The  angle 
of  attack,  a,  in  the  aerobal 1 istic  system  is  the  angle  in  the  vertical  plane 
measured  from  the  x-axis  to  the  velocity  vector;  the  angle  of  sideslip,  B  ,  • 
is  the  angle  in  the  horizontal  plane  measured  also  from J:he  x-axis  to  the 
velocity  vector.  For  the  small  yaw  angles  considered,  a  -  -n^.  and  b  *  - n y ^ . 

It  is  convenient  to  combine  a  and  s  into  a  single  complex  variable: 

|  =  8  +  i  a  -  -  (nyE  +  i  n^)  (2.1) 

The  fluid  pressure  and  viscous  forces  on  the  cavity  surfaces  produced  by 
the  motion  give  rise  to  a  moment  on  the  projectile.  The  spin-decelerating 
component,  M,  o-  ,  is  zero  to  the  approximation  considered  here;  the  other  com¬ 
ponents  can  be  represented  in  complex  form,  M|  y  +  if^T)  .  We  shall  consider 

only  the  liquid  moment  acting  on  the  projectile;  the  liqutd  moment  can  be 
added  to  the  other  moments  acting  on  shell  or  gyroscope  as  required.  The 
differential  equation  of  yawing  motion  is* 

I  d 2X/ dt 2  -  i$  Ix  dc/dt  +  I  M  X  =  i  (MLy  +  iML~).  (2.2) 

The  quantity  I  is  the  moment  of  inertia  of  the  empty  axi symmetric  shell  :-about 

A 

its  longitudinal  axis.  I  is  the  transverse  moment  of  inertia  of  the  empty 
shell  about  its  center  of  gravity.  The  spin  rate  of  the  shell  is  $,  which  is 

A 

taken  to  be  positive  in  this  work,  and  t  is  time.  The  term  I  M  X  is  an  aero¬ 
dynamic  moment  for  a  projectile.  For  a  gyroscope  this  term  is  a  gravitational 
moment  arising  from  the  separation  of  center  of  gravity  and  pivot  point,  and 
in  most  experiments  is  zero. 


*  This  is  Eq.  (2.4)  of  Reference  10  with,  only  the  liquid  moment  on  the 
right-hand  side. 
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In  general  there  is  an  interaction  between  the  motion  of  the  projectile 
and  the  liquid  motion.  Here  we  shall  specify  the  motion  of  the  projectile. 
In  particular  the  cylinder  is  nutating  with  constant  frequency  and  exponen¬ 
tially-growing  yaw: 


X  •  (Koe£T*t)e1*T*t)  =  Kje'*1  =  K0elf^, 


(2.3) 


where 


Kj  i  KoeE,,t.  =  x*t,  f  5  (1-1  e)t. 


(2.4) 


Here  KQ  is  the  magnitude  of  the  yaw  at  time  t  =  0,  t  is  the  nutational  , 

frequency  divided  by  and  e  is  a  yaw  growth  rate  or  decay  per  nutational 
cycle.  Also  is  the  yaw  amplitude,  and  <j>^  is  the  angular  orientation*  of 

the  x-axis  in  the  x,  y,  z  system  as  shown  in  Figure  1. 

The  motion  of  the  projectile  enters  the  flow  problem  via  the  boundary 
conditions.  Under  the  assumption  that  the  flow  is  in  phase  with  the  motion  of 
the  shell,  the  pressure  disturbance  will  have  the  time  dependence  of  F.q. 

(2.3),  and  consequently  the  liquid  moment  will  also  have  this  form.  A 
nondimensional  liquid  moment  coefficient,  C^M,  is  new  defined:** 

MLY  +  1  MLZ  =  \  a  *  T  CLM  Kle  l»  (2,S) 


where  mL  is  the  mass  of  the  liquid  and  a  is  the  radius  of  the  cylinder  cross- 

section.  CLM  is  a  complex  quantity  whose  real  part  represents  a  moment  that 

changes  the  yaw  angle,  and  whose  imaginary  part  represents  a  moment  that 
changes  the  nutation  rate.  Thus:  . 

CLM  E  CLSM  +  1  CLIK*  ^2*6^  ’ 


where  CLS(V,  and  CLjM  represent  the  "liquid  side  moment"  and  "liquid  in-plane 

moment,"  respectively.  As  in  Reference  9,  we  shall  concentra  te  our  attention . 
on  CLSM,  the  liquid  side  moment. 


*  For  simplicity  the  angle  of  attack  is  assumed  to  he  initially  zero  and  the 
,.uvjlc  of  sideslip  to  he  initially  positive;  i.e.,  4>  of  Reference  10  is 
zero . 

*  See  Eq.  (2.7)  in  Reference  10. 
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When  the  forcing  moment  in  Eq.  (2.5),  produced  by  the  motion  specified  in 
Eq.  (2.3),  is  inserted  into  the  yaw  equation,  Eg.  (2.2),  it  is  seen  that  the 
;  of  Eq.  (2.3)  is  a  solution  to  Eq.  (2.2)  for  a  restricted  set  of  f's,  namely 
those  satisfying  the  functional  equation 


I/2  -  Ixf  -  IyM/|2  -  -i  ( 2 upa 4C)  t  CLM  (f;  Re,  c/a). 


(2.71 


Here  c  is  the  half-height  of  the  cylinder,  p  is  the  density  of  the  liquid,  a 
Re  is  Reynolds  number  defined  by 


Re  =  a  |/v, 

where  v  is  kinematic  viscosity  of  the  liquid. 


(2.8 1 


According  to  the  Stewartson-Wedemeyer  theory4’6,  C^M  of  Eq.  (2.7)  is 
negligibly  small  except  near  resonance.  A  resonance  condition  will  generally 
occur  when  tn  Cp,  where  xn  $  is  the  nutational  frequency  of  the  empty  shell 
ana  Cp  $  is  a  natural  inertial  frequency  of  the  rotating  liquid. 


T  —  I  +  ( I 
n  x  v  x 


!y  M/I2)  1/2  (2  Iy). 


(2.9) 


For  t  *  CR,  CLM  can  be  approximated  by  the  first  term  of  the  Laurent  series  of 
a  function  with  a  simple  pole4’6 

CLM  =  D/(f  -  CR), 


where  the  residue,  0,  depends  on  the  parameters  of  the  problem. 


III.  COORDINATE  AND  VELOCITY  TRANSFORMATIONS 

4  The  flew  problem  is  stated  and  solved  in  terms  of  the  inertial  cylindrical 

coordinates  x,  r,  0  (where  y  =  r  cos  8,  z  =  r  sin  a).  However,  the  pressure 
and  viscous  forces  are  integrated  over  constant  x  and  r  surfaces  to  obtain 
moments  (where  y  =  r  cos  S,  z  =  r  sin  9).  Also,  the  original  statement  of 
boundary  conditions  is  made  in  terms  of  x,  r,  a  coordinates,  tivjs,  it  is  use¬ 
ful  to  have  the  tranf ormati ons  between  the  two  coordinate  systems.  From  Eq . 

|  (9)  of  Reference  9,  applicable  for  small  K(J , 

^  _  ?. 

r  =  r  -  K.  (x  -  xl  cos  (  j.,-0)  +  0  (K  ) 

1  1  u 

5  =  9  -  K,  [ (x  -  z)/r]  sin  (^-o)  +  0  (kq‘) 

<  2 

x  •-  x  +  r  COS  (d^-8)  +  0  (Kq  ), 

11 


I 


where  l  is  the  x  (and  x)  coordinate  of  the  pivot  point.  •  All  the  above  terms 
are  nondimensional;  lengths  and  distances  are  nondimensional ized  by  a.  The-  . 
transformation  may  also  be  expressed  as 


y  =  r  cos  e  +  (Kj  cos  ^ )  (x  -  l) 
z  =  r  sin  e  +  (K^  sin  <^)  (x  -  z) 
x  =  x  -  K^r  cos  ( ^  -  e) . 


(3.2) 


We  define  the  flow  to  be  a  small  disturbance  to  a  basic  flow,  which  is 
taken  to  be  solid-body  rotation  in  an  unyawed  cylinder.  The  Navi er-Stokes 
equations  are  linearized  to  produce  the  perturbation  equations.*  The  flow 
variables  are  the  radial,  azimuthal,  and  axial  velocity  components ,- and  pres- 
sure,  given  here  in  nondimensional  form: 


★ 

-  K0U> 


v„s  =  v  -  K0V, 


WNS  '  W 


K0«. 


-  P  -  KqP*  (3.3) 


NS  “ 

| 

0 

The  symbols  u^^,  vNS,  wNS,  and  p^<.  represent  the  total  values,  solutions  Of 

the  linearized  Navi er-Stokes  equations.  U,  V,  W,  and  P  are  the  basic  un-„ 

^ 

disturbed  variables;  and  u,  v,  w,  and  p  are  perturbation  variables**  of  order 
one.***  For  solid-body  rotation,  the  basic  flow  is 


U  =  0,  V  =  r,  W=  0,  P  =  ( 1  /2 )  r 2  +  const. 


(3.4) 


The  velocity  components  are  nondimensional  ized  by  a$>,  and  pressure  by  pa2<£2. 

In  the  aerobal 1 i Stic  system, the  radial,  azimuthal,  and  axial  velocity 

★  ★  ★ 

perturbation  components  are  denoted  by  u,  v,  and  w.  The  velocity  transfor¬ 
mation,  from  Eq.  (A. 2)  of  Reference  9,  is  given  by 


*  'These  are  Eq.  (3)  in  Reference  3. 

**  The  negative  signs  in  Eq.  (3.3)  were  employed  to  comply  with  the 
nomenclature  of  Reference  10. 

**  7 he  symbols  u,,^t  v^,q,  revlace  the  symbols  u,  v,  wy  p  in 

Eq.  (10)  of  Reference  .9. 
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( x-  £)e  9  [etcos  (<^-0)  +  (1-t)  sin  (<^-0)]  +  0(Kq) 


( x  -  i)e  9  [-(1-t)  cos  +  e  t  sin  (<^-0)]  +  0(Kq) 


(3.5) 


w  =  w  -  r  eeT<^  [-ex  cos  (<^-0)  -  (1-t)  sin  (^-0)]  +  0(Kq).  ^ 

The  tilde  superscripts  can  be  dropped  from  the  second  terms  of  the  right-hand 
sides  of  Eq.  (3.1)  and  (3.5)  without  changing  the  order  of  error. 


IV.  FLOW  SOLUTIUN  AND  WALL  FORCES 


A.  Flow  Solution 


Chapter  Ill  of  Reference  9  treats  the  flow  problem  in  detail;  we  extract 
from  it  what  is  needed  here.  The  flow  variables  are  shown  in  Eqs.  (3.3)  and 
(3.4);  the  perturbed  flow  solution  is 


u  =  Real  [u(r,x)  exp  |i(f$t-0)}] 


v  =  Real  jv(r.x)  exp  |1  (fit-0))] 
w  =  Real  [w(r,x)  exp  |i(f#t-0)}] 


(^.1) 


p  =  Real  [p ( r , x )  exp  ji(fJt-O)}] 

where  u,  v,  w,  and  g  are  complex  functions.  These  are  expressed  as  the  sums 
of  two  solutions: 


usu  +  1  u  =  u  +  u  , 

'R  '1  “p  "H 


V  =  V  +  1  V  =  V  +  V 

"R  '1  "p  “H 


(4.2) 


W  S  W  +  1  w  =  w  +  w  , 

■R  "I  "p  "h 


G  =  G  +  i  P  =  G  +  p  * 
R  I  p  "H 


where 


u  =  -i  ((l-f)2/(l  +  f)  Jx  ♦  i(l-f)i,  v  =  -t(l-f)2/(l  +  f)  ]x  +  (1-f )  £ 

"p  "p 


(4.3) 


w  =  i  (1-f  )r, 
P 


P  =  -(1-f  )2rx  +  ( 1  - f  ^ )  ir 


is  a  particular  solution.  The  sub-H  quantities  are  solutions  of  the  following 
equations,  where  subscripts  denote  partial  differentiation  and  the  sub-H  is 
omitted  for  clarity: 

ru+u-1v  +  rw=0  j 

"r  “  '  "x 


i  ( f - 1  )u  -  2v  =  -g  +  (1/Re)[u  +  u  /r  -  2u/r2  +  u  +  2i  v/ r 2  ] 

r  “rr  ”r  "  "xx 

S( 4 . 4 ) 

i (f - 1  )v  +  2u  =  i[)/r  +  (1/Re)[v  +  v  /r  -  2v / r 2  +  v  -  2 i u / r 2 ] 

"rr  "r  "  "xx 


i (f-l)w  =  -p  +  ( 1/Re) [w  +  w  /r  -  w/r 2  +  w  ]  . 

-x  "rr  "r  "  "xx 


These  quantities  further  satisfy  the  boundary  conditions  at  the  sidewall 


u  (r=l)  =  -i (2f (1-f )/( 1+f )  ]x 
"H 

v  (r=l)  =  - [2f  (1-f )/(l+f ) ]x  (4.5) 

"H  ( 

w  (r=l)  =  0, 

H 

and  at  the  endwall  (where  c  =  c/a) 


(x=c) 

-  -i  [2f ( 1-f )/( 1+f)  ;c, 

u  (x=-C)  = 

-u  (x=c) 

(4.6a) 

"H 

"H 

(x=c) 

=  -  (2f ( 1-f )/( 1+f) )c. 

v  (x=-c)  = 

:  -V  (x=c) 

(4.6b) 

"H 

"H 

(x=c) 

=  w  (x=-c)  =  0, 

(4.6c) 

"H  "H 


plus  boundary  conditions  at  r  =  0. 


Ci 


1 


•• 


e 

s 


As  explained  in  Reference  9,  a  modal  solution  (separation  of  variables)  is 
required,  but  It  cannot  satisfy  Eqs.  (4.6a)  and  (4.6b).  The  need  to  drop  two 
of  the  three  emiwall  conditions  implies  that  we  have  a  singular  perturbation 
problem;  i.a,,  we  must  insert  a  boundary-layer  or  "inner"  solution  to  sati  sfy 
these.  The  technique  of  matched  asymptotic  expansions  is  used  to  treat  this 
problem.  The  solution  u  ,  v  ,  etc.,  is  decomposed  into  an  outer  solution, 

"H  "H 

valid  away  from  the  endwall,  and  an  inner  solution,  valid  near  the  endwall. 
lhe  expansions  for  each  are  determined  to  certain  orders  in  the  small  parame- 

ter  Re  and  a  composite  solution  is  formed.  The  velocity  gradients  at  the 
endwall  needed  for  the  shear  force  are  obtained  from  the  composite  solution. 

In  Reference  10  the  viscous  correction  of  Wedemeyer,b  originally  used  to 
correct  eigenvalues,  is  applied  to  correct  velocities  and  pressure.  Although 
the  formalism  of  matched  asymptotic  expansions  is  not  used  there,  the  basic 
idea  of  obtaining  a  corrected  flow  is  carried  out.  It  is  called  the  "invis- 
cid  flow"  in  Reference  10  even  though  it  depends  on  Re;  it  would  correspond  to 
what  is  called  outer  flow  here. 


The  outer  solution,  designated  by  u,  v,  w,  p,  satisfies  Eqs.  (4.4),  (4.5), 
and  (4.6c).  The  modal  form  of  the  outer  solution  is 


u 


uk(1 


sin  x , 


(r)  sin 


v 


w  = 


(r)  cos  Akx, 


(r)  sin  A^x , 


(4.7) 


where  u^,  vk,  w^,  and  p^  are  complex  functions  of  r;  they  are  solutions  to  the 

ordinary  differential  equations,  Eqs.  (33),  with  boundary  conditions,  Eqs. 

(37)  and  (41),  of  Reference  9.  This  is  the  usual  normal  mode  solution  with 

A,  =  k ti/2c  ,  where  k  is  an  odd  integer.  In  Reference  9  it  is  shown  that  the 
outer  solution  is  determined  by  the  single  condition  of  no  flow  through  the 

1/2 

endwall,  Eq.  (4.6c),  to  the  order  0  (Re~  )•,  Eqs.  (4.6a)  and  (4.6b)  are  not 
used  or  satisfied.  The  accuracy  of  the  solution  to  this  order  Is  unaccepta¬ 
ble.  To  Improve  it,  the  second  term  of  the  outer  solution  must  be  obtained, 
which,  in  turn,  requires  the  first  term  of  the  inner  solution  and  appropriate 
matching  of  the  outer  and  inner  solutions.  As  shown  in  Reference  5, this  pro¬ 
cess  yields  the  one  boundary  condition 


w  *  6c  w^  =  0 


at  x  =  i  c 


(4.H) 
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to  the  order  0  (Re“  ).  (it  has  since  been  shown  that  this  boundary  condition 

is  correct  to  0  (Re”  ^  ).)  The  form  of  the  solution  is  still  Eq.  (4.7),  but 
the  \y  are  complex,  determined  from  the  eigenvalue  relation 


cos 


A^c  +  A^  5c  sin  A^c  =  0, 


where  5c  is  given  by  the  following  sequence: 
a  =  2”1/2Re1/2(l-i)  (3-f ) 1/2 
6=2"  ^2Re  ^2(l  +  i)  (1+f) 1,7 2 


6c  = 


7^  l  1  '  TT  )  +  h  (l  +  T2 


(4.9) 


(4.10) 


The  complex  square  roots  are  chosen  to  be  the  ones  that  make  the  real  parts  of 
a  and  0  positive.  When  6c  /c<<1,  A.  can  be  approximated  by 


\  =  (k  tt)/  [2(c  -  6c)  ].  (k  odd) 


The  first  term  of  the  inner  solution  is  determined  by  the  boundary  layer 
equations,  the  no-slip  conditions  on  the  endwall,  and  boundary  conditions  at 
the  boundary-! ayer  edge  derived  by  matching.  These  are  given  in  Reference  9, 
pages  53  and  54. 

The  theory  of  matched  asymptotic  expansions,  MAE,  is  the  proper  technique 
for  dealing  with  these  problems;  in  order  to  explain  the  results  for  shear 
force  given  in  this  report  the  discussion  of  MAE  in  Reference  9  must  be 
augmented.  Usually  MAE  are  used  to  obtain  an  analytic  solution  to  a  problem 
with,  possibly,  some  numerical  integration  required;  this  can  be  done  for  the 
present  problem,  at  least  in  principle.  Of  course,  it  becomes  increasingly 
tedious  to  obtain  higher  order  terms.  Since  the  problem  of  Eqs.  (4.4)  -  (4.6) 
is  linear,  application  of  MAE  here  is  simpler  than  for  many  other  cases  to 
which  it  is  applied.  Advantages  of  MAE  are  the  systematized  approach  to  the 
terms  in  the  expansions,  the  clear  distinction  between  inner  and  outer 
solutions  and  the  matching  of  these. 

Here  the  formalism  of  MAE  is  used  to  (1)  distinguish  the  outer  solution, 
determined  as  above,  (2)  rationalize  the  use  of  only  one  boundary  condition, 
Eq.  (4.6c),  rather  than  three  at  the  endwall,  and  (3)  derive  Eq.  (4.8).  The 
solution,  Eq.  (4.7),  with  appropriate  A^  would  be  the  exact  outer  solution 

(requiring  only  ordinary  differential  equations  to  be  integrated  numerically), 
except  that  Eq.  (4.8)  is  not  exact.  In  the  final  form  of  the  solution, the 


analytically  determined  outer  flow  is  replaced  by  the  solution  of  Eqs.  (4.7) 
and  (4.9).  This  step  leads  to  more  accurate  results  at  low  Re  though, 
strictly  speaking,  it  is  not  part  of  the  theory. 


The  first  term  of  the  inner  solution  is  given,  in  the  notation  of 
Appendix  D,  Reference  9,  by 

ui  (y,r)  =  u°  (c,  r)  +  (i/2)  [Ae*^  -  Se‘3y] 

vi  (y,r)  =  v°  (c,  r)  -  (1/2)  [Ae'0y  +  Be‘By  ] 

wi  (y,r)  =  Re  /  (av°/3x)x==  |  5c  -  y  +  [( 1+f )/2a(  1-f )  ]  exp  (-ay)  - 


1) 


[( 3-f )/20(l-f )  ]  exp  (-By)  |  A 


where  y  =  c  -  x 
and 

A(r)  -  vj  *  iu2  j 

at  x  =  c.  (4.12) 

B(r)  =  v°  -  iu°  +  4f  [( 1-f )/( 1+f )  j  c  ) 

The  functions  u°,  v°,  and  w °  are  defined  in  Reference  9,  page  53.  The  expres¬ 
sions  for  u.j  and  in  Eq.  (4.11)  would  be  the  same  as  the  boundary-layer 

solutions  of  Wedemeyer6  if  the  inviscid  terms  in  the  latter  were  evaluated  at 
the  endwall . 


From  the  inner  and  outer  expansions  a  single  expansion,  uniformly  valid  in 
the  inner  and  outer  regions,  can  be  constructed;  it  is  called  the  composite 
expansion.  Using  a  1-term  inner  and  a  2-term  outer  expansi on, the  composite 
expansions  are; 


UC  =  uo  (c-y  ,r)  +  Re“  1/2  u°q] 
VC  =  (c-y ,r )  +  Re’ 1/2  v°j 
Wr  =  w°  (c-y  ,r)  +  Re"  1/2  w° 


(c-y,r)  +  (i/2)  [Ae_oty  -  Be  'By] 
(c-y  ,r)  -  (1/2)  [Ae_oy  +  Be"Sy  j 
(c-y ,r) , 


(4.13) 


where  the  functions  u°^,  are  defined  in  Reference  9,  page  53.  The 

gradients  of  and  Vq,  (3/3y)  q  =  -  (3/3x)x_=,  are  required  to  obtain  the 
endwal  1  shear.  The  matching  process  shows  that  I~ ( 3x ) x- 5  1  of  uoi»  v°» 
v®,  are  all  zero.  Therefore,  ^ 


(  3uc/3y)y=0  =  -  (i/2)  [oA  -eB ] 
(  3yc/ 3y  )y=Q  =  (1/2)  [aA  +  83] 


(4.14) 


Thus, these  gradients  at  the  endwal 1  could  be  computed  from  just  the  inner 
solution. 


At  the  sidewall,  the  no-slip  boundary  conditions  are  satisfied  by  the  solu¬ 
tion  to  Eqs.  (4.4)  and  (4.5).  The  necessary  gradients  at  the  sidewall  are 
obtained  directly  from  the  modal  solution. 


B.  Shear  Forces 

Forces  on  the  sjrface  of  the  cylinder  are  obtained  from  the  3x3  stress 
tensor,  which  gives  the  force  on  the  fluid.  (See,  e.g.,  Reference  13,  paqe  53, 
where  1  is  axial  coordinate.)  Since  we  want  the  force  on  the  cylinder,  the 
sign  of  the  stress  tensor  is  changed  from  that  of  Reference  13.  The  elements 


of  the  tensor  are 

xrr  -  P  -  Ko  [p  -  (2/Re) 3u/ 3r ]  (4.15a) 

T00  =  P  -  Kq  [p  -  2  (u  +  3*/ 3  0)/ (Re  r)  ]  (4.15b) 

txx  =  P  -  Kq  [p  -  (2/Re)  (3w/Jx)j  (4.15c) 

Ve  =  V  =  (V**)  (r  3  (v/r)/3r  +  ' '*,l8)/r ]  (4.15d) 

Trx  =  Txr  =  <Ko/Re)  l^/3x  +  3*'/3rJ  ( 4 . 1 5e ) 

t9x  =  Txa:  (K0/Re)  (3*/3x  +  ^1/r-  3“/3e] 


( 4  - 1 5  f ) 


The  nondimensional  force  on  an  element  of  the  cylinder  surface  with  non- 
dimensional  area  dA  will  be  denoted  by 

dF  =  (dFp,  dFflf  dFx),  (4.16) 

where  dF  .  dF„,  dFv  are  the  force  components  in  the  radial  ,  azimuthal  ,  and 
axial  directions,  respectively,  in  the  earth-fixed  frame.  Force  is  non- 
dimensional  ized  by  pa4<p",  where  p  is  density  of  the  liquid,  and  area  is  non- 
dimensional  ized  by  a:.  Then 


dFr  - 

(  Trr 

q  + 
sr 

Tr  8  9  e  + 

Vx  9x>  dA 

(4.17a) 

dFe  = 

(t0T 

9r  + 

Teo  9  e  + 

T0x  9x)  dA 

(4.17b) 

dFx  = 

(Txr 

9r  + 

Tx  8  9  e  + 

Txx  9x)  dA‘ 

(4.17c) 

where  gp,  gg,  gx  are  the  radial,  azimuthal,  and  axial  components  in  the 

inertial  frame  of  a  unit  vector  (directed  outward  from  the  container;  i.e., 
away  from  the  fluid)  normal  to  the  element  of  wall  surface.  If  the  surface  is 
given  in  the  form  G(r,  9,  x)  =  const,  then 


.  3G/  Dr 

-  TvCr  ’ 


-  ( 1/r ) 3G/ 39 


(4.18) 


wnere  |  vG|  =  +  [(3G/3r)2  +  (1/r2)  (3G/3e)2+  ( 3G/3x)  2]1//2 

At  the  sidewall,  r(r,  e,  x)  =  1,  and  at  the  endwalls,  x(r,  a,  x)  =  +  c. 
From  Eqs.  (3.1)  and  (3.2)  we  obtain  unit  normals,  N  =  (g  ,  g  ,  g  )  to  sidewall 
and  endwall,  accurate  to  order  <o  :  rax 

Nside  =  (1*  K1  Ux-£)/r]  sin  K\  cos  [V0])  (4.19a) 

Ntop  =  <"K1  cos  f^i"0]*  _Ki  sin  [ ♦i - 9 ] »  1)  (4.19b) 

^bottom  '  "Ntop  (  ,*‘A) 


The  nondimensional  surface  element  areas  are 

dA-  =  d9  dx,  dA-  =  =  r  dr  da.  (4.20) 

r=l  x=±c  v  ' 

Combining  Eqs.  (4.15),  (4.17),  (4.19),  and  (4.20),  we  obtain  the  surface 
forces  correct  to  order  K  . 
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At  the  sldewal 1  (r  =  1) 


dFr/(d?  dx)  =  P  -  Kq  [p  -  (2/Re)  3u/3r 


dFe  Ko 


d  e  dx  Re 


*  ★ 

r  i_  I  + 1  iH. 

3r  r  r  38 


+  KjP  ~  sin  ( ♦1- 0) 


(4.21) 


dF  / ( d 9  dx)  =  (K  /Re)  [3u/3x  3w/3rJ  +  K.P  cos  ($,-0). 

X  0  i  1 

At  the  top  wall  (x  =  c) 

dF  /(r  dr  d  9)  =  -K^P  cos  (<^-9)  +  (K  /Re)  [3u/3x  +  3w/3rl 
dF0/(r  dr  da)  =  -K^P  sin  (^-0)  +  (K^/Re)  [3v/3x  +  (3w/30)/r] 
dFx/(r  dr"  d9)  =  P  -  K  [p  ••  (2/Re)  3w/3x]. 


(4.22) 


At  this  point  we  introduce  the  boundary-layer  assumptions,  namely,  that 
the  tangential  gradients  of  the  velocity  components  at  a  surface  are  negli¬ 
gible,  and  that  the  normal  gradient  of  the  normal  component  is  also  negligi¬ 
ble.  The  orders  of  magnitude  of  the  velocity  gradients  at  x  =  c  in  Eq.  (4.22) 
can  be  determined  explicitly  from  the  MAE  results,  Eq.  (4.13),  and  results  of 
the  matching: 

3u/3x  =  0  (Re  /2)  3w/ 3r  =  0  (Re~1//2) 

3v/3x  =  0  (Re  '  )  3w/39  =  0  (Re  '  ) 

3w/  3x  =  0  ( 1 )  . 


This  approximation  is  not  necessary  for  the  further  development  of  the 
theory;  it  is  made  here  for  convenience.  It  does  restrict  the  applicability 
of  the  moment  calculations  to  high  (as  yet  undefined)  Reynolds  numbers.  How¬ 
ever,  all  the  terms  in  the  stress  r°nsor  are  available  from  the  solution  to 
the  flow  problem;  their  contribution  to  thu  moment  may  be  significant  at  low 
Reynolds  numbers.  Actually,  a  restriction  to  high  Reynolds  numbers  has 
already  oeen  introduced  by  the  boundary  condition,  Eq.  (4.8). 
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Equations  (4.21)  and  (4.22)  reduce  to 


dF p / ( d  9  dx)  -  P  -KQp 
dF g/(d 9  dx)  =  (Kq /Re) 
dFx/(d9  dx)  =  (Ko/Re) 
at  the  si dewa 1 1 ,  and 


[r  3(v/r)/3rj  +  K^P  [_(x-*)/r]  sin  (<^-0) 

3w/3r  +  K,P  cos  ( - 9) 

1  1 


(4.23a) 


(4.23b) 


dF  /(r  dr  d9) 

dF  ,/(r  dr  d  9) 
o 

dF  /(r  dr  d9) 


-K^P  cos 

Uj-9)  +  (KQ/Re) 

3u/  3x 

(4.24a) 

-K,P  sin 

1 

Uj9)  +  (K0/Re) 

3v/3x 

(4,24b) 

★ 

P  -K  p 

A  ' 

(4.24c) 

at  tne  top  wall.  The  expressions  for  force  components  on  the  bottom  wall  are 
the  negatives  of  those  for  the  components  at  the  top  wall,  and  -c  replaces 

c. 


V.  EVALUATION  OF  LIQUID  MOMENT 
A.  Expression  for  Liquid  Moment 

We  wish  to  determine  the  moment  produced  by  the  liquid  on  the  spinning  and 
nutating  shell,  namely,  M|_y  +  i  of  Eq.  (2.5).  We  shall  evaluate  the 

moment  about  the  center  of  gravity  of  the  projectilejn  the  x,  y,  z  system. 
Details  need  be  shown  for  only  one  component,  say  since  the  form  of  Eq. 

(2.5)  indicates  that  both  M^y  and  are  determined  by  and  CLjf,  of  Eq. 

(2.6) .  In  rectangular  coordi nates, the  moment  on  an  element  of  wall  surface  is 
thej/ector  product  of  (1)  the  radius  vector  relative  to  pivot  point  (x-£, 

y,  z)  and  (2)  the  force  dF  =  (dF-,  dF-,  dF-).  The  particular  component  that 
we  treat  isM^:  x  ' 


dM 


LZ 


=  F35t2  [(x-  l)  dF-  -  y  dF-], 


(5.1) 


where  dF-  and  dF-  are  the  components  of  dF  in  the  x  and  y  directions,  respec¬ 
tively.  Unit  vectors  inj.be  x  and  y  directions  are  found  by  taking  the 
normalized  gradients  of  x  =  const,  and  y  =  const,  of  Eq.  (3.2)  in  the  manner 
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of  Eq.  (4.18);  the  scalar  products  of  these  unit  vectors  with  dF  of  Eq.  (4.16) 
yield  the  two  force  components 


dF~  =  dFr  cos  6  -  dFg  sin  6  +  dF^  K ^  cos  ^  +  0  (KQ2)  (5.2a) 

dF?  =  "dFr  K1  cos  («r0)  ■<JF0  Ki  sin  Up0)  +  dFx  +  0  V)  (5.2b) 

B.  Sidewall  Moment 

The  element  of  moment,  dM^,  in  Eq.  (5.1)  is  evaluated  at  r  =  1,  and  will 
be  denoted  by  dM^*  5y  5qs.  (3.4)  and  (3.1), 

P  (r  =  l)  =  const.  -K^  (x-£)  cos  (^-9)  i-  0  (Kq2).  (5.3) 

Application  of  Eqs.  (4.23)  and  (5.2),  yields 


dF- 

- =  (P-K  p)  cos  0  - 

dS  dx  0 


!°  ri-X,  K,P  iHil  sin 
He  3r  r  1  r 


X 

frf  -&) 


sin  6  +  0  (K  2) 
o 


f  C  1  \ 
\  o  .-4  I 


d F ~/ ( d 9  dx)  =  (K  /Re)  3w/3r  +  0  (Ko2) 


We  apply  Eqs.  (5.3)  and  (5.4)  and  integrate  Eq.  (5.1)  at  r  =  1,  noting 
that  the  constant  part  of  P  makes  no  contribution  to  the  integral,  leaving 
only  first  order  terms  in  Kq .  Thus, the  variables  r,  9,  x  may  be  replaced  by 

r,  0,  x  and  the  integral  evaluated  at  r  ■■  1  without  changing  the  first  order 
accuracy.  Integration  yields 


j  ^L/ (K0  oa  54>2)  =  -  j  j  (x-i)  |  p  *  (x-  a)  cos  (^-e)  J 


(  3w/ar)/Re  J 


co?  0  do  dx 


(5.5) 


f  (*  [”  (x-t)  f  3v/3r  -  v  \  ]  sin  e  0  3  dx  +  0  (K  2). 


(1/Re) 


Further  manipulations  lead  to  the  following  formula,  with  application  of  tqs. 
(2.51,  (2.6),  (4.1),  and  (4.2),  plus  the  boundary  condition  v(r=l)  =  -(1-f)  x 
(x-x); 


ML2-L/(2irpa4c$2T  K1 )  -  C^LSM^pL  sin  ^  +  C^LIM0L  cos  ^ 

C(LSM)VL  Sin  H  +  C ( L I M ) V L  C0S  'h 


(5.0) 


where  integrals  are  evaluated  at  r=l,  and 


C(LSK)PL  =  ^2t  j  Pj  dx 


(5.7a; 


C(LIM)PL  =  "^2  T 


if  c 

1  (x- a)  [p  +  (x-  A)  ]  hx 


(5.7b) 


'(LSM)VL 


I  ( «  (  3?  r 5:  1  ) 

7Tr7  )_=  \  ?r  *  <*-*>  L^M1-T)  (Hj 


(5.7c) 


c  --  -1-.-  r 

-IM)VL  2  t  c  Re  J_- 


c  ‘  3w  f  8v  1  ) 

7  "R  +  (x-t)  ~1  +  et  (x-£)  \ 

(  3r  L  3r  J  J 


+  E  T  (X-  *>  dx. 


( 5 . 7  d ) 


In  tne  labeling  of  the  moment  coefficients  of  Eqs.  (5.5)  and  (5.7),  the  LSM 
and  LIM  designations  are  defined  in  Chapter  II  (Eq.  (2.5)  et  seq.),  P 
indicates  pressure,  V  indicates  viscous  wall  shear,  and  the  final  L  (lateral) 
desi grates  si dewal 1 . 

C .  Endwall  Moment 

the  element  of  moment,  d’-!^*  in  Eq.  (5.1),  evaluated  at  x  r-  c  and 
x  -  -c ,  will  be  denoted  by  dM  ~T  and  dM  ~ ,  respectively.  The  element  of 

L  Z.  I  L  l.  d 

total  endwall  moment  is 

d\?E  •  ’  *Var  (b-8) 
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Application  of  Eqs.  (4.24)  and  (5.2)  yields 


5J 

« 


ts 

s> 


dF~/(r  dr  dS) 
d F ~/  (r  dr  d  6) 


(K  /Re)  [(3u/3x)  cos  6  -  (3v/3x)  sin  Q  ]  +  0  (Kq  2) 

p  -  kop  "  o  V>- 


(5.9) 


By  Eqs.  (3.4)  and  (3.1) 


?  (x  =  c)  =  fn(r)  -  K1  f  (c-i)  cos  (^-0)  +  0  (KQ2) 

P  (x  =  -c)  =  fn(r)  -  K1  r  (-c -&)  cos  (^-9)  +  0  (K02)* 


(5.10) 


We  obtain  dX,  ^  and  dM^  separately  using  Eqs.  (5.1),  (5.9),  arid  (5.10). 

The  equations  for  the  flow  solution  in  Chapter  IV  show  that  3u/'3x  and  3v/3x 
are  even  functions  of  x,  leading  to  some  cancellations  when  the  top  and  bottom 
wall  moments  are  added  to  produce 


1  ,  2  n 


M.  f 

)c  ,  0 


(2  c/Re) 


[2c  eeTfpL  r  cos  ( $^-6)  r  p  (c)  -p  ( -c )  ]  r 2  cos  e  d 9  dr  + 

(S.ll) 


1  /•  2  7T 

0  ^0 


[(3u/3x)_  COS  9  -  ( 3v/ 3x) _  Sin  0:  r  d9  dr. 
c  C 


The  r,  9,  x  may  be  replaced  by  r,  6,  x  in  Eq.  (5.11)  and  integrals  evaluated 
at  x  -  i  without  affectiny  the  first  order  approximation  in  K  . 

Moment  coefficients  analogous  to  those  for  the  sidewall  will  now  be 
defined  for  the  endwalls: 


1l2e'/'‘~  17(36  '"T21  ) 


C(LSM)PE  S1n  *1  "  ''(LIM)PE 


^  cos  4^  + 


(5.12) 


(LSM)VE  S'n  h  +  °(LIM)VE 


-  cos 


pl* 


where 
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C(LSM)PE  =  ‘  J  {  Pj  ^  “  Pj  }  p2  dr 


(5.13a) 


:(LIM)PE  M1/(2t£)!  (£/2)*J‘{pr 


(c)  -  p  (-c)  l  r2  dr 

•R  J 


(5.1  3;.^ 


'(L  SM)7 


..,r  =  -  [l/(  t  Re)  ]  Real  j  j  3(v-iu)/3x|  ,  r  dr 


/5  in,- 

y  t  A  -/  V. 


'(LIM)VE 


wp  =  -  ,l/(  t  Re)  I  mag 


I  ^  |  3(v-iu)/3xj-; 


(5.13d) 


The  total  side  moment  coefficients  due  to  pressure  and  shear  stress, 
respectively,  are 

C(ISM)P  =  C(LSM)PL  +  C(LSM)PR  *  I 


C ( L  SM ) V  =  C(LSM)VL  +  C(LSM)VE’ 


(5.i4; 


and  finally,  the  total  side  moment  coefficient  is 


CISM  =  C(LSM)P  +  C(LSM)V 


(5.15) 


The  term  '3  (v  -  iu)/3xj  _  in  Eq.  (5.13)  is  evaluated  from  the  complete 

x=c 

solution,  Eq.  (4.2),  u  =  u  +  u  and  v  =  v  +  v  .  According  to  the  method  of 

"p  *H  "  “p  ”H 

solution  described  in  Chapter  IV,  u  and  v  are  determined  by  the  composite 

"H  "H 

solutions  uc  and  v^.  Therefore,  to  the  order  of  approximation  in  Eq.  (’>.13), 


v  -  i  u  =  v  -  i  u  +  vc'’1uc 

"p  ’p 


(5. In) 


■  2x(l-f )  2/(  1  +  f )  +  2  £( 1-f )  +  v°  -  iu£  +  Re-1/2  (v^  -  iu°J  -  Be' 
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making  use  of  Eqs.  (4.3)  and  (4.13)  and  the  definition  of  Bin  Eq.  (4.10). 
The  terms  from  the  outer  solution,  Including  those  in  B,  see  Eq.  (4.12),  are 
not  evaluated  analytically;  they  are  obtained  from  Eqs.  (4.7)  and  (4.9),  as 
explained  in  Chapter  IV.  Thus  v-iu  can  be  written 


v_iu=v  -iu  +  v  -  1u  - 

“p  -P 


[(v-iu)  „  +  ^4f(l-f)/(l+f)  J  c  ]  e’®\ 


Another  useful  form  is  obtained  by  expressing  this  In  terms  of 

v  +  v  -  i (u+u  1 : 

-p  -p 


V  -  iu  *  (v+v  )  -  1(u+u  )  -  [(v+v  )  =  -  1  (u+u  ) 


X  =  c 


X  =  c 


2 ( 1  - f )  (c-A)  ]  e“6y. 


(5.17) 


(5.18) 


The  gradient  term  in  Eqs.  (5.13c)  and  (5.13d)  is  then 

[<3(v  -  iu)/ctx]  =  =|d[(v+vp)  -  i(u+up)  ]/dx  j.  =  - 
X  ~  0  X  “  c 


(5.19) 


0[(v+v  )  =  -  i(u+u  )  =  +  2 ( 1-f )  (c- ?.)] 


X  =  c 


X  =  c 


Empirically  it  is  found  that  the  first  term  on  the  right-hand  side  of  Eq. 

(5.19)  is  small  compared  to  the  second  term  and  is  neglected  when  the  gradient 

is  evalua'  ^d;  recall  that  8  =  0(Re1/,z}.  A  strict  estimate  of  the  order  of 
magnitude  of  these  two  terms  in  Eq.  (5.19)  is  not  strai ghtforward  because  we 
have  used  u  and  v,  obtained  from  Eq.  (4.7),  rather  than  the  outer  solutions. 

If  the  gradient  is  obtained  from  Eq.  (5.16),  the  result  is 


[3  (v-iu)/ ax]  =  -2(  1-f )  2/(  1+f ) 

x=c 


B  [v°  -  iu°  +  4f (1-f )/(l  +  f )]  _ 

x-c 


(5.20) 


in  which  it  is  clear  that  the  first  term  is  0(1)  and  the  second  is  0(Re'*'/,?). 
This  result  is  reflected  in  the  conclusion  about  the  terms  in  Eq.  (5.19) 
stated  above. 


VI.  NUMERICAL  RESULTS:  MOMENT  COEFFICIENT 


A.  Effect  of  Wall  Shear 


In  Reference  10  it  was  shown  that  viscous  shear  could  have  a  significant 
effect  on  the  liquid  moment  at  moderately  high  Reynolds  numbers  (-'40,000). 
The  pressure  moments  on  sidewall  and  endwalls  were  opposed  to  each  other  in 
all  computations  performed  ( C^ lSM)P  >  ®  anc^  ^(LSM)V  ^  ^ e  the  shear 

moments  on  side  and  end  walls  were  frequently  of  the  same  sign.  Thus, the 
partial  cancellation  of  the  pressure  moments  and  the  reinforcement  of  the 
shear  moments  at  times  brought  the  two  within  the  same  order  of  magnitude. 


Figures  2  and  3  illustrate  the  effect  of  viscous  shear.  At  t  =  0.19  and 
Re  =  4xl04,  e.g.,  the  magnitude  of  shear  moment  is  greater  than  25%  of  pres¬ 
sure  moment,  certainly  a  nonnegl igible  contribution.  The  difference  between 
L SM )P  and  ^LSM  3ncreases  markedly  as  Re  decreases.  The  qualitative  behavior 

of  side  moment  can  also  be  affected  by  shear.  Thus,  at  t=  0.10  in  Figure  3, 
the  pressure  moment  coefficient  increases  monotonical ly  as  Re  decreases  to 
10  \  while  the  total  moment  coefficient  peaks  at  Re^-6300  and  changes  sign  at 
Re-2500. 


In  Figure  2,  has  a  maximum  near  the  ei genf requency  of  the  2nd 

radial,  4th  axial  inertial  mode  (n=l,  k=7)  of  the  liquid  as  predicted  by 
theory.3  The  magnitude  of  shear  moment,  ! Cj_ sm )V  I  a^so  aPPears  t0  have  a 

maximum,  though  not  as  pronounced,  in  this  region. 

B.  Comparison  with  Results  of  Murphy's  Method 

In  Reference  10  Murphy  solves  the  inviscid  perturbation  equations  of 
Stewartson  together  with  the  vi scous-corrected  boundary  conditions  of 
Wedemeyer  at  sidewall  and  endwalls,  and  integrates  the  wall  forces  to  obtain 
moments  without  assuming  that  the  coning  motion  is  nearly  in  resonance  with 
inertial  oscillations  of  the  liquid.  Comparisons  between  pressure  moments 
obtained  by  Murphy's  method  and  our  method  are  given  in  Reference  9,  where  it 
is  shown  that  agreement  is  good  at  very  high  Reynolds  numbers  but  deteriorates 
with  decreasing  Reynolds  number. 

Figure  4  shows  a  comparison  of  the  total  moment  coefficients  for  the  case 
c/a  =  4.291,  e  =  0.  For  clarity  the  curves  for  Re  =  10s  and  106  have  been 
omitted;  the  curves  from  the  two  methods  are  practically  coincident  on  the 
scale  of  this  figure.  For  Re  =  10u  agreement  is  still  fairly  good;  but  for 
Re  =  103  the  discrepancies  are  large,  and  fo>  Re  =  10 2  the  two  methods  give 
Clsm's  of  opposite  sign.  There  is  no  discernable  peak  in  for  Re  = 

103  and  102  over  the  range  of  t  considered.  The  separate  pressure  and  viscous 
moment  coefficients  for  Re  -  105  and  103  are  shown  plotted  against  t  in  Figure 
5.  At  Re  =  10 5  the  pairs  of  curves  of  the  two  theories  are  practically 
coincident;  at  Re  =  10  3  and  r=  0.10,  the  ( i  SM  )P ' s  d3^er  by  30%  and  tne 

C( L SM ) V '  S  a  ^actor 
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C.  Moments  on  Sidewall  and  Endwalls 


Figures  6  and  7  illustrate  the  behavior  of  the  separate  components  of  the 
side  moment  coefficient  for  a  cylinder  of  c/a  =  3.126,  executing  motions  of 
e  =  0.02  and  t  *  0.020,  0.045.  When  Reynolds  number  is  varied  by  changing 
only  v,  as  in  the  experiments  of  Reference  14,  the  moment  varies  directly  as 
the  moment  coefficient.  It  Is  seen  then  that  pressure  moments  of  both 
Murphy's  and  present  theory  vary  nonmonotonical ly  with  Reynolds  number,  both 
on  the  sidewall  and  the  endwalls.  A  major  source  of  the  difference  between  the 
results  of  the  two  theories  lies  in  the  methods  of  obtaining  the  outer  flow, 
or,  as  designated  in  Reference  10,  the  "inviscid"  flew  (see  p.  15,  paragraph 
before  Eq.  (4.7)).  The  results  in  Figure  6  indicate  that  this  difference  is 
small;  the  location  of  the  peak  is  a  function  of  the  nutational  frequency. 

In  Figure  7  the  variations  of  shear  moments  on  the  sidewall  and  endwalls 
are  presented;  they  are  not  monotonic  functions  of  Re  and  are  negative  in  many 
instances  so  that  they  tend  to  damp  yaw.  From  the  natures  of  the  present 
theory  and  that  of  Reference  10  the  results  from  the  two  should  approach  each 
other  for  large  Re.  The  results  in  Figure  7,  and  those  for  other  values  of  t,  | 
show  that  this  occurs  for  2xl04  <  Re  <  3.2xl05  over  the  range  .02  <  t  <  .045.  f 
These  are  surprisingly  large  values  of  Re  for  an  asymptotic  approach  of  the  4 

results  of  the  two  theories.  The  differences  in  moment  on  the  endwall.  Figure  -f 
7b,  are  relatively  small.  The  gradients  there  are  obtained  in  the  same  way  in  •' 

both  theories;  however,  the  outer  flows  are  different.  The  differences  in  f 

moments  on  the  sidewall,  Figure  7a,  are  large  for  Re<104.  The  methods  of  i 

obtaining  the  gradients  on  the  sidewall  are  quite  different  in  the  two  4 

theories.  The  only  significant  difference  in  the  results  from  the  two  § 

theories  is  found  in  Figure  7a.  For  the  same  Re  and  c/a,  the  yaw  growth  rate  f 
(shown  in  Figure  11)  computed  from  the  present  theory  agrees  better  with  f 

experimental  data  than  that  computed  from  Reference-  10. 

The  variation  of  shear  moment  with  Re  on  either  the  sidewall  or  endwall  is  4 

_  l  /  2  ■* 

not  a  simple  power  law.  On  the  endwall,  the  shear  moment  would  vary  as  Re 
if  the  analytically  determined  outer  flow  were  used  to  obtain  the  gradient  in  T 
Eq.  (5.19),  as  explained  in  the  discussion  on  Eqs.  (5.19)  and  (5.20).  The  ,s 
nonmonotonic  variation  shown  in  Figure  7b  precludes  such  a  variation;  a  power  I 
law  variation  would  plot  as  negative  exponential  .in  Figure  7b.  The  same  con-  - 
elusion  holds  for  the  sidewall,  Figure  7a,  except,  that  for  t  =  0.020  (and 
0.030,  not  shown)  the  shear  moment  is  monotonic  with  Re.  However,  if  a  power  \ 

law  fit,  Re”n,  to  the  sidewall  shear  moment  is  tried,  0.55  <  n  <  0.81  is  : 

obtained  for  5xl03  <  Re  <  5x10s.  Note  that  in  Reference  10,  p.  34,  it  is 

.1/2  l 

stated  that  all  shear  moment  coefficients  vary  as  Re  .  The  results  of 
Figure  7  shew  that  this  cannot  be  the  case. 


1-1.  W. !’.  D’Amico  and  M.C.  Miller ,  "Flight  Instability  Produced  by  a  Rapidly 
: '.pinning.  Highly  Viscous  Liquid.,  "Journal  of  Spacecraft  and ■  Rockets , 

Vnl.  ltd  January -February  1979,  pp".  G2-F4 . 
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VII.  YAW  GROWTH  RATES:  COMPARISON  OF  THEORY  ANO  EXPERIMENT 


Ultimately,  the  validity  of  the  theory  must  be  gauged  by  comparison  of 
theoretical  and  experimental  results.  Measurements  of  nutational  frequencies 
and  yaw  growth  rates  can  be  made  In  gyroscope  experiments.  The  apparatus, 
operations,  and  accuracy  of  these  experiments  are  discussed  in  References  11 
and  12;  they  are  thus  far  the  only  type  of  experiments  available  for  compari¬ 
son  with  the  theory.  The  gyroscope  can  be  used  because  the  equation 
describing  its  oscillatory  motion  is  analogous  to  that  describing  the  yawing 
motion  of  a  projectile.  The  relationship  between  gyroscope  and  projectile 
motions  is  discussed  in  Chapter  2  of  Reference  2.  The  moments  of  inertia 

I  and  I  in  Eqs.  (2.2)  and  (2.7)  necessarily  include  parts  of  the  appa- 

*  y  *  - 

ratus.  The  term  1^  H  (  in  Eq.  (2.2)  is  a  gravitational  moment  arising  from 

the  separation  of  center  of  mass  and  pivot  point.  It  will  be  zero  here  since 
these  two  positions  are  essentially  coincident;  thus  will  be  zero.  The 
theoretical  quantities,  a  and  <_ ,  are  obtained  by  solving  Eq.  (2.7)  for 
f  s  (l-ie)-r.  The  imaginary  part  of  Eq.  (2.7)  can  be  written 

etn  [1  +  {(T"Tn)/Tn)}j  L1  +{(T-Tn)/T}  (2*pa-c)  CL$M/ly.  (7.1) 

It  was  found  that  |  t- -tn  |  /  t  0.06  for  the  experimental  cases;  in  fact, 

I  T"  Tnl  /  r  t  •(^3  in  most  instances.  Thus  e  varied  roughly  as  (C^/I^). 

In  the  experiments  circular  cylinders  were  filled  with  liquid.  Table  1 
(page  30)  shows  the  seven  cases  for  which  measurements  were  taken.  The  first 
two  cases  are  the  ones  shown  in  Figure  8  of  Reference  12:  Ix  was  varied  by 

adding  flat  metal  rings  around  the  cylinder,  and  corresponding  changes  in  Iy 

were  made  by  moving  a  counterweight  to  adjust  the  position  of  the  center  of 
gravity  of  the  gyroscope.  Neither  1  nor  tn  was  recorded;  it  was  necessary  to 

estimate  I  by  a  process  which  Involved  application  of  the  Stewartson- 
Wedemeyer  theory  (p.  19  of  Reference  12).  In  the  last  five  cases,  I  was  kept 
constant  and  I  was  varied  by  moving  a  weight  along  a  threaded  shaft  coin¬ 
ciding  with  the  longitudinal  axis  of  the  cylinder11.  The  values  of  I  and 
I  used  in  the  runs  were  chosen  so  that  the  empty-shell  nutational  frequen¬ 
cies,  (  -  I x / 1 y )  ,1 ay  a  ran9e  values  covering  the  eigenfrequencies 
CR  shown  in  the  table.  Reynolds  number  was  varied  by  changing  tht  liquids, 

which  were  silicon  oils  of  differing  viscosities.  The  estimated  nominal 
values  of  I  are  4.14xl9fe  g  cm2  for  Case  1  and  4.33xl06  g  cm2  for  Case  2. 

The  n  and  k  identify  the  inertial  mode  whose  nondimensional  frequency  is 
CR  (-frequency/.j.) . 

In  Table  1  the  nominal  c/a  is  that  value  quoted  in  the  appropriate 
reference,11*1-'-  In  Reference  10  the  term  "fitted"  c/a  was  introduced.  As 
used  here,  fitted  c/a  is  determined  in  the  following  way.  Uncertainties  in  the 
measurement  of  cylinder  dimensions  give  rise  to  experimental  errors  in  c/a 


TA8LE  1.  LIST  OF  EXPERIMENTAL  CASES 


# 

Re 

IBM 

Nomi nal 
c/a 

n, 

■ 

Ref. 

1 

5.20-105 

3.164 

1, 

3 

12 

2 

9.0OxlO3 

3.152 

1, 

3 

12 

3 

5.21xl03 

3. 140 

3.126 

1, 

3 

★ 

4 

. 

l.OlxlO1 

3.130 

3.126 

1, 

3 

★ 

5 

1.24x10“ 

1.047 

1.042 

1, 

1 

11 

6 

2.4°  '10 3 

1.047 

1.042 

1, 

1 

11 

7 

1.26x10'' 

1.047 

1.042 

1, 

1 

* 

» 

(Fitted  c/a) 
CR 

p(g/cm‘:) 

I  (g  cm2) 

A 

a  (cm) 

1 

0.0486 

0.818 

Variabl  e 

3.153 

2 

0.0515 

0.960 

Variable 

3.153 

3 

0.0495 

0.972 

8.23*105 

4.121 

4 

0.0532 

0.972 

1.08xl06 

4.121 

5 

0.0437 

0.966 

7.94x10= 

6.359 

6 

0.0425 

0.972 

i.05xl0c' 

6.359 

7 

0.0415 

0.974 

1.05xl0£ 

6.359 

which  may  be  as  large  as  ±1%.*  Values  of  c/a  lying  within  the  error  bound 
about  the  nominal  value  can  therefore  be  used  in  the  calculation.  The  fitted 
value  is  that  which  gives  the  best  agreement  with  the  data,  judged  sub¬ 
jectively.  It  will  be  seen  that  the  ei  vs  t  relationship  is  extremely 
sensitive  to  c/a. 

Results  are  presented  in  Figures  8-14  in  the  form  of  plots  of  yaw  growth 
rate  versus  nutational  frequency.  In  ail  cases  theoretical  curves  are  drawn 
for  fitted  values  of  c/a,  and  in  two  cases  (Figures  10  and  12)  also  for  the 

nominal  c/a.  It  is  evident  from  the  seven  cases  that  agreement  between  theory 

and  experiment  is  better  for  the  c/a  =  1.0  than  for  the  c/a  =  3.1  cases.  It 

is  not  understood  why  this  is  so;  it  may  be  related  to  the  fact  that  resonance 

is  excited  for  a  simpler  axial  mode  (k=l)  for  c/a  -  1.0  than  for  c/a  =  3.1, 
where  the  k=3  mode  resonates.** 

In  Case  1,  Figure  8,  the  two  theories  agree  with  each  other  at  the  high 
Reynolds  number,  as  expected.  However,  the  agreement  with  experiment  is  not 
as  good  as  anticipated.  The  uncertainty  in  I  previously  mentioned  is  partly 

responsible.  Comparison  of  peak  locations  is  not  possible  because  of  lack  of 
sufficient  experimental  points  to  delineate  the  maximum  clearly  and  because 
of  scatter  in  the  data.  In  Case  2,  Figure  9,  discrepancies  between  theory  and 
experiment  are  small,  although  percentage  errors  are  large.  In  Case  3,  Figure 
10,  the  agreement  is  poor  for  the  lower  range  of  nutation  rates.  Overall, 
present  theory  results  show  better  agreement  with  experimental  data  for  Cases 
2,  3,  and  4  (Figures  9,  10,  11)  than  do  results  of  the  theory  of  Reference  10; 
it  is  expected  that  the  two  theories  would  disagree  at  lower  Re. 

We  consider  the  three  cases  for  the  c/a  -  1.0  cylinder.  For  Re  =  12,400, 
Figure  12,  the  data  and  the  two  theories  agree  quite  well.  For  Re  =  2400, 
Figure  13,  the  present  theory  agrees  somewhat  better  with  experiment  than  does 
the  theory  of  Reference  10  in  the  prediction  of  x  [U*)max]  =  xm;  for  Re  = 

1260,  Figure  14,  it  gives  the  better  overall  fit  to  the  measurements.  The 
results  in  Figures  12-14  for  Re  =  12,400,  2400,  and  1260  show  that  xm  for 

nominal  c/a  exceeds  CR  by  21%,  42%,  and  58%  (referred  to  x  ) ,  respectively,  a 

monotonic  increase  with  decreasing  Re.  These  differences  are  all  greater  than 
the  corresponding  differences  for  the  c/a  =3.1  results.  The  magnitude  of 
( xm-CR)  provides  a  measure  of  the  departure  of  the  Stewartson-lJedemeyer 

appro  i max ' on  from  our  theory.  The  present  theory  and  that  of  Reference  10, 
which  is  an  improvement  on  the  Stewartson-Wedemeyer  theory,  yield  essentially 
the  same  results  for  the  c/a  =1.0  cases.  Evidently,  both  c/a  and  Re  determine 
the  differences  in  the  results  of  the  two  theories. 


*  W.r.  D'Amico,  private  aorrsrunicatioK . 

**  The  mode  r.yjnber,  k,  occurs  in  the  flow  solution;  see  Tog.  (4.7)  und  (i.j). 


The  calculated  yaw  growth  rate  Is  quite  sensitive  to  small  changes  In 
c/a.  This  can  be  seen  In  Figure  10,  Case  3.  The  nominal  value  of  c/a  gives  a 
Tm  which  Is  13%  less  than  that  for  the  data,  whereas  the  fitted  value  gives  a 

tm  3%  less.  The  calculated  (e-r)max  are  31%  and  11%  greater  than  the  experi¬ 
mental  maximum  for  the  nominal  and  fitted  values  of  c/a,  respectively.  For 
3.126  <  c/a  <  3.140,  the  calculated  r  and  (£t)_.v  would  lie  between  the  two 

111  mu  /S 

curves  shown  in  Figure  10.  Thus,  no  choice  of  c/a  will  give  a  clearly  superior 
result. 

Sensitivity  of  the  er  vs  r  curves  to  changes  in  aspect  ratio  is  further 
illustrated  in  Table  2  for  the  parameters  of  Case  3.  The  results  show  that  a 
0.8%  change  in  c/a  produces  a  32%  change  in  (et)max  and  a  13%  change  in  t  . 

Even  over  this  small  range  of  c/a, both  quantities  depart  noticeably  from  a 
linear  variation.  This  sensitivity  has  important  implications  for  the  theory 
vis-a-vis  results  from  laboratory  experiments  and  field  firings.  With  some 
care  it  is  feasible  to  control  c/a  to  within  tolerances  of,  say,  0.5%,  or  even 
less,  in  laboratory  experiments.  Deformations  of  the  cylinder  would  have  to 
be  accounted  for;  some  possible  causes  are  clamping  of  cylinder,  compression 
of  the  liquid  upon  installation  of  the  top,  and  temperature  changes.  In  the 
experiments  quoted  here,  it  appears  that  c/a  was  known  to  an  accuracy  of 
+0.5%.  It  is  probably  not  feasible  to  control  c/a  to  that  accuracy  in  field 
firings,  '’ypical  manufacturing  processes  allow  the  internal  dimensions  of  the 
cylinder  to  vary  by  much  more  than  that.  The  application  of  the  theory  to 
such  cases  is  then  questionable. 


TABLE  2.  MAGNITUDE  AND  LOCATION  OF  (et)max 


c/a 

Tm  E  T  E(e*)Max] 

3.126 

1.358  *  10'4 

0.046 

3.140 

1.600  x  10‘4 

0.050 

3.150 

1.788  x  10‘4 

0.052 

VIII.  CONCLUSIONS 

We  have  developed  a  method  of  computing  the  moment  exerted  by  the  spun-up 
liquid  on  the  casing  of  a  filled  shell  that  is  spinning  and  nutating.  In 
addition,  we  are  able  to  predict  the  nutational  and  yaw  growth  rates  of  the 
projectile's  angular  motion. 

The  output  yielos,  separately  and  in  combination,  pressure  and  viscous 
shear  moments  on  sidewall  and  endwalls.  The  applicability  of  the  method  is 
restricted  to  small  angles  of  attack  because  of  linearization  of  the  Navier- 
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Stokes  equations  and  to  late  times  in  the  flight  history  because  of  the  as¬ 
sumption  of  initially  unperturbed  solid-body  rotation.  There  is,  in  addition, 
a  Reynolds  number  limitation  resulting  from  the  presence  of  an  error  of 
.1 

order  0  (Re  ).  The  permissible  smallness  of  Reynolds  number  has  not  been 
definitely  determined;  there  is  an  indication  that  it  depends  on  the  oscilla¬ 
tion  mode  primarily  being  excited. 

Limited  parameter  studies  of  moment  coefficients  indicated  that  sidewall 
and  endwall  pressure  moments  opposed  each  other  and  that  viscous  shear  moment 
was  often  not  negligible.  Otherwise,  no  simple  trends  were  discerned  relating 
the  relative  contributions  (sidewall,  endwall,  pressure,  shear)  to  the  side 
moment  with  the  various  parameters  of  the  problem  (Re,  c/a,  t,  e). 

Comparison  of  side  moment  coefficients  with  those  of  Murphy10  showed  good 
agreement  for  high  Reynolds  numbers,  but  increasing  divergence  with  decreasing 
Reynolds  number.  For  the  sidewall  viscous  shear  moment  the  relative  discrep¬ 
ancies  were  large  for  Re  <  1,000  (Figure  7a).  However,  the  corresponding 
discrepancies  in  yaw  growth  rate  were  not  necessarily  large  also.  For 
example,  the  parameters  of  Case  1  (Figure  8),  substituted  into  Eq.  (7.1),  led 
to  the  following  relation  between  the  discrepancies  a  (ex  )  and  a  (CICm): 

4  Ui  n)  .0.0012  MCLSH).  n  LSH 

Yaw  growth  rate  outputs  from  this  theory  were  compared  with  results  of 
Murphy's  theory  and  with  measurements  from  gyroscope  experiments.  At  high  Re 
the  differences  between  the  two  theories  were  small  but  increased  with  de¬ 
creasing  Re.  Further  experiments  covering  the  parameter  space  in  more  detail 
are  needed  to  provide  a  better  assessment  of  the  theory.  All  the  experiments 
treated  here  have  Re  <  12,400  except  for  Case  1,  which  has  an  uncertainty 
in  I  .  Thus,  the  most  obvious  need  is  for  experimental  results  for  Re  > 

50,000.  In  projectile  firings  values  of  Re, up  to  several  million  are  common¬ 
place. 

The  sensitivity  of  yaw  growth  rate  to  small  changes  in  c/a,  discussed  at 
the  end  of  Section  VII,  is  a  significant  result.  This  effect  has  been  known 
for  some  time;  the  theories  of  Reference  10  and  the  present  paper  have  pro¬ 
vided  a  definitive  demonstration  of  it. 
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igure  1.  Diagrams  of  Coordinates  and  Cylinder. 
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Figure  2. 


Pressure  Moment  Coefficients  (Sidewall,  Endwall ,  Total) 
Total  Viscous  Shear  Moment  Coefficient,  and  Total 
Moment  Coefficient  vs  Mutational  Frequency  (Re=40,000, 
c/a=4 .291 ,  s=0.0,  i=0). 
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Figure  3, 


log10  Re 

Variation  of  Pressure  Side-Moment  Coefficient  and  Total 
Side-Moment  Coefficient  with  Reynolds  Number  for  Fixed 
Nutational  Frequency  (c/a=4.291,  c=0.0,  v=0). 


Frequency:  Comparison  of  Results  of 
Reference  10  Method  and  Present  Method 
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Pressure  and  Viscous  Shear  Side-Moment  Coefficient 
vs  Nutational  Frequency  for  Re  =  10'  and  1 0 b : 
Comparison  of  Results  of  Reference  10  Method  and 
Present  Method  (c/a=4.291,  c=0.0,  i-O) . 
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Figure  7.  Viscous  Shear  Side-Moment  Coefficient  on  (a)  Sidewall 
and  (b)  Endwalls:  Comparison  of  Results  of  Reference 
10  Method  and  Present  Method  (c/a=3.126,  c=0.O2,  >.=0) 
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Figure  8.  Yaw  Growth  Rate  vs  Nutational  Frequency,  Case  1: 

Comparison  of  Experimental  Results  with  Those  of 
Reference  10  Method  and  Present  Method  (Re  = 
5.20x10s,  c/a=3. 154 — Fitted) . 


Figure  9.  Yaw  Growth  Rate  vs  Nutational  Frequency,  Case  2: 

Comparison  of  Experimental  Results  with  Those  of 
Reference  10  Method  and  Present  Method  (Re  = 
9.00xl03,  c/a=3.152--Fitted) . 
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Figure  10.  Yaw  Growth  Rate  vs  Nutational  Frequency,  Case  3: 

Comparison  of  Experimental  Results  with  Those 
of  Reference  10  Method  and  Present  Method 
(Re  =  5.21x103,  c/a=3. 140— Fitted,  c/a=3 . 126-- 
Nominal ) . 


T 

Figure  11.  Yaw  Growth  Rate  vs  Nutational  Frequency,  Case  4: 

Comparison  of  Experimental  Results  with  Those  of 
Reference  10  Method  and  Present  Method  (Re  = 
1.01x10’,  c/a=3.130--Fittea) . 


T 

Figure  12.  Yaw  Growth  Rate  vs  Nutational  Frequency,  Case  5: 

Comparison  of  Experimental  Results  with  Those  of 
Reference  10  Method  and  Present  Method 
(Re  =  1.24xlO\  c/a=l. 047— Fitted,  c/a=1.042  — 
Nominal ) . 
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Figure  13.  Yaw  Growth  Rate  vs  Nutational  Frequency,  Case  6: 

Comparison  of  Experimental  Results  with  Those  of 
Reference  10  Method  and  Present  Method 
(Re  =  2.40xl03,  c/a=l. 047— Fitted). 
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14.  Yaw  Growth  Rate  vs  Nutational  Frequency,  Case  7 
Comparison  of  Experimental  Results  with  Those 
Reference  10  Method  and  Present  Method 
(Re  =  1.26xl03,  c/a=l .047--Fitted) . 
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CR 

CLM 

CLIM 

CLSM 
C(L IM)P 

C(LSM)P 
C  (L IM)  V 
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C(LIM)PE 

C(LSM)PE 

C(LIM)PL 

C(LSM)PL 

C(LIM)VE 

C(LSM)VE 

C(LIM)VL 

C(LSM)VL 
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cross-sectional  radius  of  cylinder  [cm] 
functions  of  r  defined  in  Eq.  (4.12) 
half-height  of  cylinder  [cm] 

=  c/a ,  aspect  ratio 

natural  oscillation  frequency  of  rotating  liquid/^ 
liquid  moment  coefficient  =  C|^M  +  i  CLjM,  Eq.  (2.5) 

liquid  in-plane  moment  coefficient  =  C^i^p  +  ^(lIM)V’ 

liquid  side  moment  coefficient,  Eqs.  (2.6)  and  (5.15) 
pressure  in-plane  moment  coefficient  =  C^jM^pL  +  Cm  ,M^p£ 

pressure  side  moment  coefficient,  Eq.  (5.14) 

viscous  shear  in-plane  moment  coefficient  = 

C(LIM)VL  +  C(LIM)VE 

viscous  shear  side  moment  coefficient,  Eq.  (5.14) 

endwall  pressure  in-plane  moment  coefficient,  Eq.  (5.13b) 

endwall  pressure  side  moment  coefficient,  Eq.  (5.13a) 

sidewall  pressure  in-plane  moment  coefficient,  Eq.  (5.7b) 

sidewall  pressure  side  moment  coefficient,  Eq.  (5.7a) 

endwall  viscous  shear  in-plane  moment  coefficient, 

Eq.  (5.13d) 

endwall  viscous  shear  side  moment  coefficient,  Eq.  (5.13c) 

sidewall  viscous  shear  in-plane  moment  coefficient, 

Eq.  (5.7d) 

sidewall  viscous  shear  side  moment  coefficient,  Eq.  (5.7c) 
nondimensional  wall  surface  area  element,  Eq.  (4.20) 
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dF 

dFr,  dF9,  dFx 

dF~,  dF~ 
y  ’  x 


nondimensional  stress  force  exerted  by  liquid  on  dA 

4,2 

[Force/(pa  $  )] 

radial,  azimuthal,  and  axial  components,  respectively,  of 
dF  [Force/(pa4$>2)] 

components  of  dF  in  y  and  x  directions,  respectively 


=  (l-ie)x,  complex  representation  of  angular  motion, 

Eq.  (2.4) 

moment  of  inertia  of  empty  shell  about  its  longitudinal 
axis  [g  cm2] 

transverse  moment  of  inertia  of  empty  shell  about  its 
center  of  gravity  [g  cm2] 

index  of  axial  eigenfunction  and  eigenvalue,  Eqs.  (4.7) 
and  (4.9) 

yaw  anolitude  at  time  t  =  0 

• 

=  KqC81^,  yaw  amplitude  at  time  t,  Eqs.  (2.3)  and  (2.4) 
nondimensional  x  (and  x)  coordinate  of  pivot  point 
mass  of  liquid  in  cylinder  =  2npa2c  [g] 


\y  ’  Ml_Z 


M  M 
LZB’  LZT 


MLZE’  mlzl 


n 

nYE  ’  nZE 


aerodynamic  (or  gravity)  moment  parameter,  Eq.  (2.2) 

y  and  z  components,  respectively,  of  liquid  moment 
[g  cm  z/s 2] 

bottom  and  top  wall  contributions,  respectively,  to 
[g  cr,.2/s2] 

endwall  and  sidewall  contributions,  respectively,  to 
[g  cm2/s 2] 

Index  of  radial  mode  for  eigenf requency ,  CR 

components  in  the  y,  z  plane  of  a  unit  vector  lying  on  the 
x-axis 


N 


si  de 


unit  vector  normal  to  sidewall;  similar  definition  for 

vP> i4-1") 
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*  z  * 


p  =  p 

-  -R 


+  1 


P 

-I 


nondimensional  r,  x  variation  of  perturbation  pressure, 
Eq.  (4.1)  [pressure/(K0pa2$2)] 


P 

-H 


PNS 


P 

-P 

* 

P 


Pk(r) 


series  solution  contribution  to  p  in  Eq.  (4.4) 

total  pressure/(pa 2$2) ,  Eq.  (3.3),  solution  to  Navier- 
Stokes  equations 

particular  solution  contribution  to  p,  Eq.  (4.3) 
perturbation  pressure/(KQpa2£2) ,  Eqs.  (3.3)  and  (4.1) 
coefficient  of  sin  in  p  series,  Eq.  (4.7) 


P 

r 

r 

Re 

t 


unperturbed  pressure/(pa2$2) ,  Eqs.  (3.3)  and  (3.4) 
radial  coordinate  in  inertial  system/a 
radial  coordinate  in  nonrotating  aerobal 1 i sti c  system/a 
Reynolds  number  =  a2$/v 
time  [s] 


u,  v,  w,  p  nondimensional  r,  x  variation  of  outer  solution  of  flow 

problem 

u,  v,  w  nondimensional  r,  x  variation  of  perturbation  velocity 

components,  Eq.  (4.1)  [velocity/(KQia) ] 


nondimensional  perturbation  velocity  components  in 
inertial  system,  Eqs.  (3.3)  and  (4.1)  [vel ocity / (KQ|a) ] 


nondimensional  radial,  azimuthal,  and  axial  perturbation 
velocity  components  in  aeroballistic  system,  Eq.  (3.5) 
[velocity/(Ko$a)] 

uf,  vr,  wr,  pr  composite  solution  contribution  to  r,  x  variation  inflow 

problem,  Eq.  (5.16) 


u  ,  v  ,  w  series  solution  contribution  to  u,  v,  w,  Eqs.  (4.2) 

"H  "H  "H  and  (4.4)  '  “  " 

u-j.  vi ,  w-,  first  term  of  inner  solution,  Eq.  (05)  of  Reference  9 


/I  7 


A  A  A 

u.  ,  v.  ,  w  coefficients  of  sin  X-x  and  cos  x.x  In  u,  v,  w  series, 

*  K  k  Eq.  (4,7)  X  k 


UNS*  VNS ’  WNS 

nondlmenslonal  radial,  azimuthal,  and  axial  velocity 
components,  respectively,  in  inertial  system  (Eq. 
(3.3)),  solution  to  linearized  Navler-Stokes 
equations  [veloc1ty/($a)] 

uo  ’  vo ’  wo  ’  Po 

first  term  of  expansion  of  outer  solution  in  Eq.  (D3) 
of  Reference  9 

u°  ,  vO.  w° 

0  O  0 

first  term  of  expansion  of  uQ ,  vQ ,  wQ  in  Eq.  (D7)  of 
Reference  9 

uol’  vol’  wol’  pol 

second  term  of  expansion  of  outer  solution  in  Eq.  (03) 
of  Reference  9 

0  0  0 

U  .  ,  V  .  ,  w  , 
ol ’  ol  ol 

first  term  of  expansion  of  uQj,  vQp  in  Eq.  (08) 

of  Reference  9 

U  ,  v  ,  w 

-p  -p  -p 

particular  solution  contribution  to  u,  v,  w,  Eqs. 

(4.2)  and  (4.3)  '  "  * 

-R  *  -I 

real  and  imaginary  parts,  respectively,  of  u 
(similar  definitions  for  vD ,  vT  and  wD ,  w‘ 

Eq.  (4.2))  "K  'K  _I 

u,  v,  w 

nondimensional  radial,  azimuthal,  and  axial  velocity 
components  of  unperturbed  flow,  Eq.  (3.4)  [velocity/ 
(♦»)] 

x,  y,  i 

nondimensional  rectangular  coordinates  in  inertial 
system  (x-axis  along  trajectory)  [length/a] 

X,  y,  z 

nondimensional  rectangular  coordinates  in 

aerobal 1 i stic  system  (x-axis  along  cylinder  axis) 
[length/a] 

y 

=  c  -  x  in  Chapter  IV  (Eq.  (4.11)  et  seq.) 

a,  6 

functions  of  f  and  Re  defined  in  Eq.  (4.10) 

a 

angle  in  vertical  plane  measured  from  the  x-axis  to 
the  velocity  vector 

8 

angle  in  horizonta1  plane  measured  from  the  x-axis  to 
the  velocity  vector 

5  C 

correction  term  in  endwall  boundary  condition,  Eq. 
(4.10) 
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=  (1/t)  x  yaw  growth  per  radian  of  nutation 

polar  angles  (azimuthal  coordinates)  in  inertial  and  aero- 
ballistic  systems,  respectively 

eigenvalue  in  the  axial  problem,  Eqs.  (4.7)  and  (4.9) 

kinematic  viscosity  of  liquid  [cm2/s] 

vector  describing  angular  motion  of  cylinder,  Fq.  (2.1) 

density  of  liquid  [g/cm3] 

nutational  frequency  of  cylinder/I 

nutational  frequency  of  empty  shell /$,  Eq.  (2.9) 

nondimensional  components  of  stress  tensor,  tq.  (4.15) 
[stress/(  pa2|2)] 

=  x$t ,  angular  orientation  of  x-axis  in  the  x,  y,  z  system 
spin  rate  of  cylinder  [rad/s],  taken  to  be  positive 
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